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The spinning-top Be star Achernar from VLTI-VINCI 
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Abstract. We report here the first observations of a rapidly rotating Be star, a Eridani, using Earth-rotation 
synthesis on the Very Large Telescope (VLT) Interferometer. Our measures correspond to a 2a/2b = 1.56 ± 0.05 
apparent oblate star, 2a and 2b being the equivalent uniform disc angular diameters in the equatorial and polar 
direction. Considering the presence of a circumstellar envelope (CSE) we argue that our measurement corresponds 
to a truly distorted star since a Eridani exhibited negligible Ha emission during the interferometric observations. 
In this framework we conclude that the commonly adopted Roche approximation (uniform rotation and centrally 
condensed mass) should not apply to a Eridani. This result opens new perspectives to basic astrophysical problems, 
such as rotationally enhanced mass loss and internal angular momentum distribution. In addition to its intimate 
relation with magnetism and pulsation, rapid rotation thus provides a key to the Be phenomenon: one of the 
outstanding non-resolved problems in stellar physics. 
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1. Introduction 

The southern star Achernar (a Eridani, HD 10144, spec- 
tral type B3Vpe) is the brightest Be star in the sky. A Be 
star is defined as a non-supergiant B type star that has 
presented episodic Balmer lines in emission (Jaschek et al. 
1981), whose origin is attributed to a CSE ejected by the 
star itself. Physical mechanisms like non-radial pulsations, 
magnetic activity, or binarity are invoked to explain the 
CSE formation of Be stars in conjunction with their funda- 
mental property of rapid rotation. Theoretically, rotation 
has several consequences on the star structure (Cassinelli 
1987). The most obvious is the geometrical deformation 
that results in a larger radius at the equator than at the 
poles. Another well established effect, known as gravity 
darkening or the von Zeipel effect for hot stars (von Zeipel 
1924), is that both surface gravity and emitted flux de- 
crease from the poles to the equator. Although well stud- 
ied in the literature, the effects of rotation have rarely 
been tested against accurate enough observations (Reiners 
& Schmitt 2003, van Belle et al. 2001), a gap bridged by 
our interferometric observations of Achernar. 
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2. Observations and data processing 

Dedicated observations of Achernar have been carried 
out during the ESO period 70, from 11 September to 
12 November 2002, with quasi-uniform time coverage, on 
the VLT Interferometer (VLTI, Glindemann et al. 2003) 
equipped with the VINCI beam combiner (Kervella et al. 
2003) . This instrument recombines the light from two tele- 
scopes in the astronomical K band, which is centered at 
2.2 /jm and covers 0.4 /im. The observable measured by 
VINCI is the squared coherence factor /i 2 of the star light. 
It is derived from the raw interferograms after photomet- 
ric calibration using a wavelet based method (Segransan et 
al. 1999). The reduction procedure is detailed by Kervella 
et al. (2003) and has successfully been applied to dwarf 
stars observations with the VLTI (Segransan et al. 2003). 
The instrumental value of n 2 is then calibrated through 
the observation of stable stars with known angular diame- 
ters. The calibrators chosen for Achernar are presented in 
Table 1. The final product of the processing is the squared 
visibility V 2 of the object for each baseline projected on 
the sky (£? pro j). V 2 is directly related to the Fourier trans- 
form of the brightness distribution of the object via the 
Zcrnike-Van Cittert theorem. For these observations, two 
interferometric baselines were used, 66 m (E0-G1; azimuth 
147°, counted from North to East) and 140 m (B3-M0; 
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Fig. 1. VLTI ground baselines for Achernar observations 
and their corresponding projections onto the sky at dif- 
ferent observing times. Left: Aerial view of VLTI ground 
baselines for the two pairs of 35 ccm siderostats used for 
Achernar observations. Color magenta represents the 66 m 
(E0-G1; azimuth 147°, counted from North to East) and 
green the 140 m (B3-M0; 58°). Right: Corresponding base- 
line projections onto the sky (-B pro j) as seen from the star. 
Note the very efficient Earth-rotation synthesis resulting 
in a nearly complete coverage in azimuth angles. 

58°), equipped with 40 cm siderostats (Fig. 1 left). Their 
orientations are almost perpendicular to each other giv- 
ing an excellent configuration for the detection of stel- 
lar asymmetries. Moreover, Earth-rotation has produced 
an efficient baseline synthesis effect (Fig. 1 right). A to- 
tal of more than 20000 interferograms were recorded on 
Achernar, and approximately as many on its calibrators, 
corresponding to more than 20 hours of integration. From 
these data, we obtained 60 individual V 2 estimates, at an 
effective wavelength of A ff = 2.175 ± 0.003 /im. 

Table 1. Relevant parameters of the calibrators of 
Achernar. 0ud is the equivalent uniform disc angular di- 
ameter. The value of 0ud for 8 Phe and \ Phe is based on 
spectrophotometry (Cohen et al. 1999), while that for a 
PsA was measured separately with the VLTI and should 
appear in a forthcoming publication from one of us (E.F.). 



Name 


Spec.type 


A c ff 


Baseline 


0UD 






(H 


(m) 


(mas) 


S Phe 


KOIIIb 


2.181 


140 


2.18 ±0.02 


a PsA 


A3V 


2.177 


140 


2.20 ±0.07 


X Phe 


K5III 


2.182 


66 


2.69 ±0.03 



3. Results 

Determination of the shape of Achernar from our set of V 2 
is not a straightforward task so that some prior assump- 



tions need to be made in order to construct an initial solu- 
tion for our observations. A convenient first approximation 
is to derive from each V 2 an equivalent uniform disc (UD) 
angular diameter 0tjd from the relation V 2 — \2J i(z) / z\ 2 . 
Here, z — 7t0tjd(oO Bproj (a) A^, Ji is the Bessel function 
of the first kind and of first order, and a is the azimuth 
angle of i? pro j at different observing times due to Earth- 
rotation. The application of this simple procedure reveals 
the extremely oblate shape of Achernar from the distribu- 
tion of 0tjd(o ! ) on an ellipse (Fig. 2). Since a, B plo j(a), 
and A c ff are known much better than 1%, the measured 
errors in V 2 are associated only to the uncertainties in 
0UD- We performed a non-linear regression fit using the 
equation of an ellipse in polar coordinates. Although this 
equation can be linearized in Cartesian coordinates, such a 
procedure was preferred to preserve the original, and sup- 
posedly Gaussian, residuals distribution as well as to cor- 
rectly determine the parameters and their expected errors. 
We find a major axis 2a — 2.53 ± 0.06 milliarcsec (mas), 
a minor axis 2b = 1.62 ± 0.01 mas, and a minor-axis ori- 
entation ao = 39° ± 1°. Note that the corresponding ratio 
2a/2b = 1.56 ±0.05 determines the equivalent star oblate- 
ness only in a first-order UD approximation. To interpret 
our data in terms of physical parameters of Achernar, a 
consistent scenario must be tailored from its basic known 
properties, so that we can safely establish the conditions 
where a coherent model can be built and discussed. 

4. Discussion 

Achernar's pronounced apparent asymmetry obtained in 
this first approximation, together with the fact that it is 
a Be star, raises the question of whether we observe the 
stellar photosphere with or without an additional contri- 
bution from a CSE. 

For example, a flattened envelope in the equatorial 
plane would increase the apparent oblateness of the star 
if it were to introduce a significant infrared (IR) excess 
with respect to the photospheric continuum. Theoretical 
models (Poeckert & Marlborough 1978) predict a rather 
low CSE contribution in the K band especially for a star 
tilted at higher inclinations, which should be our case as 
discussed below. Indeed, Yudin (2001) reported a near IR 
excess (difference between observed and standard color in- 
dices in visible and L band centered at 3.6 /im) to be E(V- 
L) = 0™2, with the same level of uncertainty. Moreover, 
this author reports a zero intrinsic polarization (p*). These 
values are significantly smaller than mean values for Be 
stars earlier than B3 ( E(V — L) > 0'. n 5 and pZ > 0.6%), 
meaning that the Achernar's CSE is weaker than in other 
known Be stars. Further, an intermediate angular diam- 
eter of our elliptical fit (Fig. 2) is compatible with the 
0ud = 1-85 ± 0.07 mas measured by Hanbury-Brown 
(1974) in the visible, in contrast to what is expected if 
the envelope were to contribute to our present IR obser- 
vations. Finally, Chauville et al. (2001) report no emission 
in the H7 line. Since the emission in H7 and in the contin- 
uum are both formed roughly in the same layer of the CSE 
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Fig. 2. Fit of an ellipse over the observed squared visi- 
bilities V 2 translated to equivalent uniform disc angular 
diameters. Each V 2 is plotted together with its symmet- 
rical value in azimuth. Magenta points are for the 66 m 
baseline and green points are for the 140 m baseline. The 
fitted ellipse results in major axis 2a = 2.53±0.06 milliarc- 
sec, minor axis 2b — 1.62 ± 0.01 milliarcsec, and minor 
axis orientation ao = 39°±1° (from North to East). The 
points distribution reveals an extremely oblate shape with 
a ratio 2a/2b= 1.56 ±0.05. 

(Ballereau et al. 1995), the contribution from the nearby 
environment of the star should be considered below the 
limit of detection. 

Of course, being classified as a Be star, Achernar can 
enhance the strength of its CSE due to episodic mass ejec- 
tions, which are generally witnessed by increased Balmer 
line emission (e.g. de Freitas Pacheco 1982). This pos- 
sibility was checked against a Ha spectrum (Leister & 
Janot-Pacheco 2003) taken in October 2002, during our 
VLTI- VINCI campaign presenting a photospheric absorp- 
tion profile. To be sure that we observed a quiescent 
Achernar we synthesized a Ha profile from our model 
(Domiciano de Souza et al. 2002) compared to the ob- 
served line. We estimated the emission to be at most 3% 
across the whole line. Such an upper limit would corre- 
spond to a CSE emitting at most 12% of the photospheric 
continuum flux, due to free-free and free-bound emission 
(Poeckert & Marlborough 1978). 

Thus, we assume hereafter that the observed asymme- 
try of Achernar mainly reflects its true photospheric dis- 
tortion with a negligible CSE contribution. Under this as- 
sumption, and using the Hipparcos distance (d = 44.1 ±1.1 
pc; Perryman et al. 1997), we derive an equatorial ra- 
dius R cq — 12.0 ± 0.4 Rq and a maximum polar radius 
R™* = 7.7±0.2 R Q , respectively from 2a and 2b obtained 



from the elliptical fit on 0tjd(o:). From simple geometrical 
considerations the actual polar radius i? po i will be smaller 
than i?™ ^ x for polar inclinations i < 90°, while R cq is 
independent of i. 

Table 2. Fundamental parameters for two limit solu- 
tion models of Achernar. From the fixed parameters and 
in the Roche approximation, the minimum polar incli- 
nation is i m ; n = 46° where V cq = V cr it — 311 km/s 
and R cq — 1.5i? po i- In addition to these parameters we 
adopted a linear limb darkening parameter from Claret 
(2000) compatible with the variable effective temperature 
and gravity over the stellar surface. The equatorial effec- 
tive temperature is much lower than the polar one due to 
the von Zeipel effect. 



Fixed parameters 


Adopted value 


Comments 


Tpoi 


20 000 K 


~ B3V star 


Mass 


6.07 M 


Harmanec (1988) 


Kq sin i 


225 km/s 


Slettebak (1982) 


Reel 


12.0 R 


this work 


Model dependent 


Values for 


Values for 


parameters 


Model A 


Model B 




9 500 K 


14 800 K 


i 


50° 


90° 


Virit 


304 km/s 


285 km/s 


V eq 


0.96V cri t 


0.79V cri t 


-Rpoi 


8.3 R 


9.5 R 



Based on these conclusions we applied our 
intcrfcrometry-oriented code (Domiciano de Souza 
et al. 2002) to Achernar. This code includes radiation 
transfer, the von Zeipel law (T c g cx g^g 5 , T e s and g e s 
being the effective temperature and gravity, respectively), 
and the Roche approximation (e.g. Roche 1837, Kopal 
1987). In this approximation and noting that the stellar 
rotation must be kept smaller than its critical value, the 
adopted projected equatorial velocity V eq sini = 225 km/s 
(Slettebak 1982) implies that i > 46°. At this limit T cff 
and g c ff both attain zero at the equator, and the surface 
equipotential first derivatives become discontinuous. 
Therefore we chose to explore a parameter space between 
the representative limit solution models A (i — 50°) 
and B (i — 90°). Table 1 summarizes the corresponding 
sets of fundamental parameters. Fig. 3 clearly shows 
that the solutions enclosed between the models A and B 
cannot reproduce the observed highly oblate ellipse. We 
also checked, with negative result, whether the situation 
would improve significantly by varying the fundamental 
parameters of Achernar in a physically reasonable range 
(mass ±1M Q , T pol ± 2000K, V cq sini ± 25 km/s). 

Thus, in absence of Ha emission making a CSE con- 
tribution unlikely to reproduce the observed oblatcncss, 
the classical assumption of Roche approximation becomes 
questionable. Deviations from this gravitational potential 
and the presence of differential rotation, both intimately 
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Fig. 3. Comparison of ratios of squared visibility curves 
between the polar and equatorial directions Vp ol /V^ q . The 
black solid curve corresponds to V^ ol /V^ q for the ellipti- 
cal fit on Achernar 's data, together with the correspond- 
ing uncertainties. Simply speaking Vp ol /V^ q somehow re- 
flects R cq /R™f, since interferometry is sensitive to the 
Fourier transform of the stellar brightness distribution. 
The colored region represents our attempt attain the black 
curve with our model for Achernar within the physically 
reasonable solutions A (orange; upper limit) and B (red; 
lower limit). This failure to reproduce the observations is 
a strong and direct indication that uniform rotation does 
not apply to rapidly rotating stars. 



related to the internal angular momentum distribution, 
should be investigated. Indeed, several differential rotation 
theories predict surface deformations stronger than that of 
uniform rotation by considering that the angular velocity 
increases towards the stellar center. Two interesting exam- 
ples are "shellular" rotation (Zahn 1992) and laws with 
angular momentum constant on cylinders (Bodcnhcimcr 
1971). In this context our result on Achernar's surface 
distortion should also impact other internal mechanisms 
like meridional circulation, turbulence, transport and dif- 
fusion of the chemical elements and angular momentum, 
increase of mass loss with rotation as well as anisotropies 
in the mass ejection and wind density from rotating stars 
(Maeder 1999, Maeder & Meynet 2000). 

Finally, the highly distorted shape of Achernar poses 
the question of Be stars rotation rate. As argued by several 
authors (e.g. Cassinclli 1987, Owocki 2003) the formation 
of out-flowing discs from Be stars remains their central 
puzzle, where rapid rotation is the crucial piece. Struve's 
(1931) original vision of a critically rotating star, eject- 
ing material from its equator, has been discarded in the 



past by observing that Be stars rotate at most 70% or 
80% of their critical velocity (typically ~ 500 km/s for a 
B0V star). However, this statistically observed limit might 
be biased by the fact that close to or beyond such veloci- 
ties the diagnosis of Doppler-broadened spectral lines fails 
to determine the rotation value due to gravity darkening 
(Owocki 2003, Townsend 1997). We believe that only di- 
rect measures of Be star photospheres by interferometry 
can overcome the challenge to prove whether these objects 
rotate close, to a few percent, of their critical velocity or 
not. This have a profound impact on dynamical models 
of Be stars CSE formation from rapid rotation combined 
to mechanisms like pulsation, radiation pressure of photo- 
spheric hot spots, or expelled plasma by magnetic flares. 
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